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Abstract. Parton cascade is a microscopic transport approach for the study of the space-time evolution 
of the Quark-Gluon Plasma produced in relativistic heavy ion collisions and its experimental manifesta- 
tions. In the following, parton cascade calculations on elliptic flow and thermalization will be discussed. 
Dynamical evolution is shown to be important for the production of elliptic flow including the scaling and 
the breaking of the scaling of elliptic flow. The degree of thermalization is estimated using both an elastic 
parton cascade and a radiative transport model. A longitudinal to transverse pressure ratio, Pl/Pt ~ 0.8, 
is shown to be expected in the central cell in central collisions. This provides information on viscous 
corrections to the ideal hydrodynamical approach. 

PACS. 24.10.Lx Monte Carlo simulations - 25.75.-q Relativistic heavy-ion collisions 
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1 Introduction 

Quantum Chromodynamics (QCD) is the theory of strong 
interactions. To understand strong interaction phenom- 
ena, people invented many QCD motivated models. One 
such model is the parton cascade model. It is a natural 
continuation of the hadron transport approach l,2[[3l[4]. 
One major advantage of this microscopic transport de- 
scription of relativistic heavy ion collisions is that it does 
not rely on the assumption of local thermal equilibrium. 
In other words, it can be used to study the equilibra- 
tion process. The parton cascade concept was introduced 
by Klaus Kinder-Geiger and Berndt Muller 0[6]. Klaus 
named his parton cascade program VNI ( "Vincent Le Cu- 
Curullo Con GiGinello" for "that little guy who plays with 
quarks and gluons" ) [7] . Many developments of VNI and 
other parton cascades 8.9. 10 . 1TUT2] have lead to new in- 
sights into Quark-Gluon Plasma production. 

Parton cascade describes the evolution of the partonic 
system by solving the Boltzmann equation. Some simpli- 
fications are introduced to make the problem tractable. 
Interactions between particles and the color field are not 
included. Particles' color degrees of freedom are not fol- 
lowed. Under these conditions, the Boltzmann equation 
can be written as 



d p d 
dt + E"dx 



f(x,p,t) = 
S(x,p, t) + C 22 + C 23 + C 32 + 



In the above equation, f(x,p, t) is the phase-space dis- 
tribution. Its time evolution depends on the source term 
S(x,p,t) and the collision terms C 22 , C 2 3, C32, etc. The 
source term describes particle production from processes 



other than direct collisions, e.g., particle production from 
strong color field or Glasma 13]. Each collision term, C mn , 
describes collisions with m incoming particles and n out- 
going particles. It is an integral of dimensions 3 x (m + n — 
1) — 4 with m phase-space distributions. In general, the 
evolution of the phase-space distribution can not be solved 
analytically. One way of solving the Boltzmann equation 
numerically is to discretize the phase space distribution. 
It can be written as a sum of contributions from point 
particles, 



f(x,p,t) 



5^(p- Pi ). 



Where particle i with weight Wi and four- momentum (pj, Ei) 
propagates in a straight line from production position Xi 
at production time t%. Then Monte Carlo method is used 
to evaluate the collision terms. 

In the following, recent results from parton cascade 
and microscopic transport will be reported. The main fo- 
cus will be on elliptic flow (Sec. [2]) and thermalization 
(Sec. [3|). The former is important for the understanding 
of experimental data and the latter offers insight into the 
macroscopic description of heavy ion collisions. The dis- 
cussion will end with a summary and outlook. 



2 Elliptic flow 

Elliptic flow measures transverse momentum anisotropy 
relative to the reaction plane. It can be characterized by 
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the second Fourier coefficient of the particle azimuthal dis- 
tribution. If a Quark-Gluon Plasma is produced, large el- 
liptic flow is expected based on the fact that partons form 
early than hadrons and the partonic equation of state is 
harder than the hadronic equation of state. Since the ellip- 
tic flow is produced early during the evolution, it reflects 
interactions of partons inside the Quark-Gluon Plasma. 
This was indeed demonstrated by Zhang, Gyulassy and 
Ko [14] . As RHIC data came out [15] , a study using more 
realistic diffuse transverse gluon geometry instead of over- 
lapping cylinders was carried out by Molnar and Gyulassy 
[TO] , It was found out that large elastic cross sections 
(on the order of 45mb) are needed in order to describe 
RHIC elliptic flow data. This certainly implies that pro- 
cesses other than perturbative elastic gluon scatterings 
can be important. Further calculations were carried out 
by Xu and Greiner [12] . They implemented the two glu- 
ons to three gluons and its inverse reaction. The Landau- 
Pomcranchuk-Migdal (LPM) effect is approximated by cut- 
ting off radiated gluons that can not form before the next 
collision. This leads to large angle radiative scatterings 
which facilitate thermalization and the buildup of elliptic 
flow. Under the assumption of local parton-hadron dual- 
ity, Xu, Greiner and Stocker showed that the calculated 
elliptic flow matches experimental data if the strong in- 
teraction coupling constant a s = 0.6 is used |16j . 

Within the elastic parton cascade framework, recent 
works have been done to investigate the relation between 
the macroscopic description and the microscopic descrip- 
tion of heavy ion collisions [TTlfTB"] . Molnar and Huovinen 
showed that the microscopic and macroscopic descriptions 
agree in the description of transverse momentum p t dif- 
ferential flow [18] . This is true not only for the case with 
a very large constant rescattering cross section, but also 
for the case with a time dependent cross section tuned 
such that the shear viscosity to entropy density ratio i]/s 
is approximately 1/ (47r). In other words, when the macro- 
scopic description, i.e., viscous hydro applies, the macro- 
scopic description is equivalent to the microscopic descrip- 
tion in describing experimental observables. When radia- 
tive processes are included, rj/s is shown to be very small. 
In particular, if a s = 0.6 is used, rj/s is close to l/(47r) 
[16 . This is what motivated many recent works on parton 
transport at minimal viscosity. One recent progress in this 
area is the demonstration of elliptic flow scaling [19] by 
Ferini, Colonna, Di Toro and Greco. It is shown that when 
freeze-out happens at an energy density value of e = 0.2 
GeV/fm 3 , the pt differential elliptic flow scales with both 
the initial spatial anisotropy e x and with the integrated 
elliptic flow {V2). However, when the system freezes out at 
e = 0.5 GeV/fm 3 , the V2/e x scaling is broken, while the 
V2/{k{v2)) still holds. This interesting result revised the 
common belief that the integrated vi is a good measure of 
the initial spatial anisotropy. The breaking of V2/e x scal- 
ing is larger for larger impact parameter. It is consistent 
with the breaking of vi/t x scaling observed experimen- 
tally [20 . This breaking gives valuable information about 
the evolution and freeze-out of the partonic system. 



3 Thermalization 

Ideal hydrodynamics is very successful in describing RHIC 
data ; 21, 22,23, 24J. One assumption of ideal hydrodynam- 
ics is local thermal equilibrium. Hydro dynamical equa- 
tions can also be used when there is local isotropy [35] . In 
this case, additional entropy production equation needs to 
be used in place of entropy conservation. Quantum Me- 
chanics requires non-zero viscosity [26 . The local equi- 
librium assumption of ideal hydrodynamics implies zero 
viscosity. Hence ideal hydrodynamics is only an effective 
description. The degree of thermalization is an important 
aspect of improving the ideal hydrodynamical description 
of relativistic heavy ion collisions [27 ] l28 l l29] . 

The evolution of bulk properties of the hot and dense 
matter produced in relativistic heavy ion collisions was 
studied recently in the frame work of the AMPT (A Multi- 
Phase Transport) model [30]. The AMPT model is a hy- 
brid model [Tll[3T][32l[33]. It uses HIJING (Heavy Ion 
Jet-INteraction Generator) [33] as the initial condition. 
The publicly available version has two options: the de- 
fault model and the string melting model. In the parton 
stage, the default model has only mini-jet gluons while 
the string melting model melts hadrons in HIJING into 
quarks and anti-quarks according to their valence struc- 
tures. The parton evolution is handled with the ZPC 
(Zhang's Parton Cascade) parton cascade model |9j. The 
default model hadronizes via the Lund string fragmen- 
tation [35] while the string melting model uses a simple 
coalescence model [361137] for hadronization. The ART (A 
Relativistic Transport) model 2,38 is used for the hadron 
transport processes. The default model is good at describ- 
ing particle distributions 31,32] and the string melting 
model is needed for the description of elliptic flow [3S] and 
HBT radii [40] . Other observables have also been studied 
using the AMPT model ^I4^l4l]l44ll431l4^l471l4^l49] . 

The bulk properties of matter in the central cell in 
central collisions are studied. In particular, the equation 
of state is characterized by the pressure over energy den- 
sity ratio P/e as a function of energy density. In the cen- 
tral cell, the equations of state of both the default model 
and the string melting model interpolate between the hard 
partonic phase and the soft hadronic phase. However, they 
differ in details, especially in the intermediate energy den- 
sity range. It is important to notice that unlike for the 
equilibrium case, the equation of state is not the whole 
story. To gauge the degree of equilibration, additional in- 
formation is needed. For central heavy ion collisions, be- 
cause of cylindrical symmetry, the longitudinal pressure 
to transverse pressure ratio, Pl/Pt, can be used to char- 
acterize the degree of thermalization. More precisely, it 
is a measure of the degree of isotropization. Fig. Q] shows 
Pl/Pt from the AMPT model. Due to the increase of 
parton production in the initial stage, the string melting 
model has a faster increase in the early stage. As the par- 
ton cross section increases, Pl/Pt goes closer to 1. The 
pressure anisotropy approaches 1, then it increases and 
crosses 1. In this particular case, Pl/Pt = 1 does not im- 
ply thermalization. It is caused by the onset of transverse 
expansion. This also shows up in the change of slope in the 
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Fig. 1. Evolution of P l /Pt from the AMPT model. 

energy density evolution. The exponent of energy density 
evolution changes as transverse expansion sets in. This 
certainly demonstrates that only partial thermalization is 
achieved. 

To study more carefully the effects of radiative pro- 
cesses, we will look at whether pressure isotropy can be 
maintained by gluon elastic and inelastic processes. For 
this purpose, we start with an initial local thermal gluon 
distribution. Gluons are produced at a proper time r = 0.5 
fm/c. The initial space-time rapidity distribution is uni- 
form with a space-time rapidity density of 1000 and the 
gluons are distributed between space-time rapidity rj s = 
—5 and +5. In the transverse direction, at the formation 
proper time, they are distributed uniformly within a ra- 
dius of 5 fm. We will start with an initial temperature 
To = 0.5 GeV. The evolution in r is studied with fixed 
grid in r] s and expanding grid in the transverse direction. 
16 test particles per real particle are used in the calcula- 
tions. Fig. Ogives the time evolution of Pl/Pt- One can 
use the free streaming curve as a reference for the effect of 
expansion. In the free streaming case, if the average lon- 
gitudinal momentum squared over the average transverse 
momentum squared, < pf > / < >, is used to measure 
the isotropy, the evolution follows exactly 1/r 2 . For the 
Pl I Pt evolution, free streaming can be approximated by 
1/t 2 , but more precisely, it evolves slightly slower than 
1/r 2 . Now we turn on 2 to 2 (two incoming gluons and 
two outgoing gluons), 2 to 3 and its inverse reaction. The 
cross sections are chosen to be isotropic to maximize equi- 
libration. C23 is fixed to be 20% of 0-2,2 on the same order 
as measured in Ref. |12j . I32 is determined by detailed 
balance. (T22 is taken to be the screened Coulomb cross 
section with the screening mass determined dynamically. 
Two values of the strong interaction coupling constant 
a s are used. The figure shows that when interactions are 
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Fig. 2. Pl/Pt evolution starting with a local thermal initial 
condition. 



turned on, Pt /Pt evolution deviates from the free stream- 
ing curve. The longitudinal to transverse pressure ratio 
becomes larger as a s increases. Due to the competition of 
expansion and thermalization, there is a minimum in the 
Pl I Pt evolution. The larger the strong coupling constant 
is, the stronger the thermalization is, and the earlier the 
time that the minimum occurs. 

It is interesting to see what happens if the parameters 
of the system change. It turns out there is a scaling law for 
the case with elastic collisions with dynamical screening. 
More precisely, the Pt / Pt evolution depends only on the 
combination a s To. This is illustrated in Fig. [3J The two 
curves in Fig. [3] have different initial thermal conditions. 
The solid one has To = 0.5 GeV and the dashed one has 
T = 1 GeV. However, they have the same a s T = 0.3 
GeV. The Pl/Pt curves agree with each other. This scal- 
ing is caused by the same initial density (no) and the same 
initial binary cross section (cr 2 2 °c olsTq/uq). This combi- 
nation leads to the same evolution of density and pressure 
isotropy. When 2 to 3 and 3 to 2 are turned on, different 
systems with different initial temperatures will evolve to- 
ward different chemical equilibrium densities. This breaks 
the a s To scaling. Fig. 2] shows the evolution of Pl/Pt 
when 2 to 3 and 3 to 2 are included. It is interesting to 
notice that though the a s T scaling is broken, the evolu- 
tion has an approximate a s scaling. This certainly demon- 
strates the important of particle number changing pro- 
cesses in thermalization. The case with a s = 0.6 gives a 
value of Pl/Pt that is around 0.8. It is consistent with 
the AMPT results with parton rescattering cross section 
a = 10 mb as shown in Fig. Q] Calculations from other 
models [50 , 51 , 52J also give comparable anisotropy. This 
level of anisotropy should be expected in the improved 
viscous hydrodynamical studies. 
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Fig. 3. a s To scaling in Pl/Pt evolution for elastic scattering 
with dynamical screening mass. 



A recent study by Fochler, Xu and Greiner [33] shows 
that parton transport may be able to give a unified de- 
scription of both jet quenching and the large pt differ- 
ential elliptic flow. This is not easy by straightforward 
application of jet quenching formalism as the expansion 
of the system makes it difficult to generate enough ellip- 
tic flow [54], Straightforward application of jet quench- 
ing formalism also faces problem with explaining both 
the light quark jet quenching and the heavy quark jet 
quenching. Several studies have shown the importance of 
transport processes, including approaches with heavy res- 
onances , the Langevin model approach [57] , effec- 
tive large elastic cross sections for heavy quarks [35] ; early 
formation and dissociation of heavy mesons [58] , and the 
enhancement of A c /D ratio in heavy ion collisions |59j . 
Partonic transport can certainly contribute to the under- 
standing of this phenomenon by including these mecha- 
nisms into the space-time evolution of the system. There 
are also recent parton transport studies on the ridge and 
Mach cone phenomena [60ll61[[62] that show the impor- 
tance of strong interactions. Further studies with partonic 
transport will help to quantitatively understand these and 
many other phenomena. 
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Fig. 4. Approximate ct s scaling in Pl/Pt evolution. 



4 Summary and outlook 

The parton cascade/microscopic transport method has been 
used to study relativistic heavy ion collisions. It has made 
important contributions to the understanding of the physics 
behind elliptic flow and thermalization. In particular, trans- 
port theory and viscous hydrodynamics agree for the de- 
scription of dense systems. Parton freezeout is seen to be 
important for elliptic flow scaling. The degree of kinetic 
equilibration depends on detailed transport processes. Ra- 
diative processes can have unique contributions to ther- 
malization. 
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